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ABSTRAC 



In this study ore three-day and two five-day Numerical 
Operational Global Atmospheric Prediction System (NCGAES) 
500mb forecasts were spectrally verified. The wavenumber 
compcnents cf the observed and forecast waves were grouped 
into planetary (wavenumbers 1-3) , long (wavenumbers 5-7) and 
medium (wavenumbers 8-12) wavenumbers for verification. The 
observed ana forecast trough-ridge (Hcvmcller) diagrams cf 
longitude versus time for each group were analysed for each 
case . 

In all three cases the most serious errors occurred in 
the planetary waves where the model often forecast errone- 
ously large cr small amplitudes. The long waves were the 
most accurately forecast group both in amplitude and phase 
speed. The medium wave amplitudes were forecast too weak 
and phase speeds were consistently too fast for this grou 
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INTROD UCTION 



Accurate and reliable numerical weather prediction is a 
major goal fcr civilian and military atmospheric studies. 
Verification studies are an important aspect cf numerical 
model development as they isolate model strengths and weak- 
nesses and identify errors which must be eliminated tc 
extend the present limits cf forecast skill. 

This thesis will investigate the Naval Operational 
Global Atmospheric Prediction System (NOGAPS) model’s abil- 
ity tc forecast the position and movement of 500mb waves in 
the mid latitudes of the Northern Hemisphere. NOGAPS, which 
became operational in September, 1982 at the Fleet Numerical 
Oceanography Center (FNOC) has b<=en a project of the Naval 
Environmental Prediction Research Facility (NEPRF) since 
1976. The icng term goal cf NOGAPS is to achieve a medium 
range (7- 1C days) numerical forecast capability fcr the U.S. 
Navy (Rosmcnd , 198 1) . Verification cf N0GAP3 is very impor- 
tant at this stace due tc the research type background cf 
the medal. The NOGAPS forecast modal is a modified form of 
the CCLA general circulation model. Consequently, it has 

not had the benefit of many prior intensive verification 
studies as have most second and third generation operational 
models. 
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Recent studies have raised serious questions concerning 
the relative accuracy of the planetary-scale (zonal wavenum- 
bers 1-2) wave motion forecasts at 500mb produced by numeri- 
cal weather prediction models. Lambert and Merilees (1978) 
concluded the planetary-scale motion predictions of a spec- 
tral numerical model were poorer with respect to persistence 
than long (wavenumbers 3-5) or medium-scale (wavenumbers 
6-10) motions. Miyakcda et al. (1972) examined the GFDL 

model performance over a series of winter forecasts and 
found similar error characteristics in the planetary-scale 
forecasts. This study will verify NOGAPS planetary wave 
forecasts and ether wave regimes, large and medium-scale. 

500mb data from both the forecast and the analysis are 
spectrally decomposed into its basic components of amplitude 
and phase than organized into planetary, long and medium- 
scale waves for separate analysis of their atmospheric 
behavioral characteristics. If one assumes spatial coupling 
(the interacting of different scales) , then model weaknesses 
in fcreca sting 'long* wave motion will be reflected in the 
accuracy of 'short' wave motion forecasts. The purpose cf 
this study is tc perform a thorough analysis cf the fore- 
casts of each wave grouping. Somerville (1980) examined the 
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planetary -seals wave forecasts cf a primitive equation 
numerical model in be th a global and hemrspheric configura- 
tion. He concluded that the global version was markedly 
more skillful than the hemispheric version especially in the 
latter part cf the five-day forecast period. While the 
present study will consider only a glcoal version of NGGASS 
the results will lend themselves to an application of Somer- 
ville’s work to NCGAPS. 

Errors in phase speeds, amplitudes, development of barc- 
clir.io systems, damping cf smaller scale features and plane- 
tary wave structure become clearer when the forecasted 
variables are analyzed over time via a Hovmoller or trough- 
ridge plot. Baumhefner and Downey (1978) utilized this 
technique to interccmpare the forecasting skills of three 
numerical weather prediction models. They constructed lon- 
gitude-time plots along a fixed latitude to better under- 
stand the time evolution cf the forecast fields. Forecasts 
which originally appeared quite good were revealed to have 
amplification errors in many of the transient systems. A 
spectral Ireakdcwn cf the Hovmcller plot in planetary 
(wavenumbers 1-3), large (wavenumbers 5-7) ar.d medium 
scale-waves (wavenumbers 8-12) were utilized in their study. 
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On? has only to consider the complexity of a typical 500mb 
wave pattern tc appreciate the benefits of this approach. 

The specific objectives cf this thesis are: 

(1) E>e velopment of a spectral verification program using 
wavenumber grouping and Hcvmoller plots for tropospheric 



height data, 

(2) Application of the program to interesting NOGAPS fore- 
cast situations at 500mb, 

(3) Use of wavenumber analysis to verify NOGAPS forecasts 
of planetary, long and medium-scale waves. 

Chapter 2 reviews the main features of the NOGAPS model. 
The pattern decomposition program is described and demon- 
strated m chapter 3. Chapter 4 presents the basic charac- 
teristics cf the case study periods and their subsequent 
atmospheric behavior. finally, chapter 5 • summarizes the 
results of the forecast comparison and makes suggestions for 
further research. 
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II 



• ncgaps summary 

T he Naval Operational Global Atmospheric Prediction Sys- 
tem (SOGAP3) was developed by the Naval Environmental Pre- 
diction Research Facility (NEPRF) to be run at the Fleet 
Numerical Cceanccraphy Center (FNOC) . The objectives of the 
NOGAPS project are to improve tropical and Southern Hemis- 
peric forecasting and achieve mid-range (5 - 10 days) global 
numerical forecast capability to support U.S. Naval opera- 
tions (Rosmcnd , 1 981 ) . 

The objective analysis method used by NOGAPS is a varia- 
tion of the successive corrections technique {Barnes, 1964) 
which discri jinates among observing systems according to 
their accuracy. While the final analysts is r.ct as sophis- 
ticated as the optimal interpolation method it requires far 
less computer resources. 

Giooal mass and wind fields are analyzed independently 
with optional data bcgusing. Cata rejection criteria are 
intentionally loose since the initialization scheme 
(3arker , 1 981) removes the gravity waves excited by bad data. 
This techniaue prevents inadvertent rejection of the good 
data associated with rapidly developing systems. 
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result; in subsequent first guess errors sufficiently large 
to cause rejection of good observational data. 

The NQGAPS model is a modified version of -he UCLA gen- 
eral crrcuiaticn model (GCM) described Dy Arakawa and Lamb 
(1977). The 2.4 by 3.0 degree horizontal grid uses scheme C 
staggering cf the variables and conserves energy and enstro- 
phy when the flew is tond iver gen t . All dependent variables 

except vertical velocity are defined in the middle of each 
of the six vertical sigma levels (Fig. 1) . 

Moisture, racist static energy and momentum are 
assumed to be well-mixed in the P3L. The equations govern- 
ing the diabetic processes are based on the method cf Dear- 
derff (1972). These equations represent the conservation cf 
mass in the PEL, temperature and specific humidity inversion 
strength and the jump in momentum across the inversion. A 
special feature of the model is the variable depth of the 
?BL ir. the lowest layer which has a capping inversion cf a 
porous material surface. This allows entrained mass ■‘■o 
deepen the PEL and removed mass (by cumulus mass flux) tc 
decrease the PEL. The fluxes at the top cf the ?BL are net 
a constant fraction of the surface fluxes but rather a func- 
tion cf the entrainment. Dissipating stratus causes a 
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Figure 1 . 



T ha Vertical Differencing of the NOGAPS Model 
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particularly large flux across the inversion for stratus 
ca pped FBI's. 

a significant modification to the GCM is the cor. strain - 
lent cf the FEI tc remain in the bottom sigma layer. This 
eliminates serious computational problems yet has net cre- 
ated any noticeable adverse effects cn model forecasts. 

NCGAPS uses the Araka wa-S chu bert (1974) cumulus parame- 
terization described by Lord (1978). Environmental air 
below the cloud top is entrained while air at the cloud top 
is detrained. Aracawa noted that large scale forcing gener- 
ates conditional stability whereas cumulus corvee- icn 
destroys it. This concept, called 'quasi-equilibrium', was 
combined with the 'cloud work' concept (Arakawa and Schu- 
bert, 1974) to produce a cloud base mass flux unique to each 
cloud type. 

The GCh radiation parameterization requires only fifteen 
percent of the model's running time. The radiation transfer 
equation (Kat ay ana , 1 974 and Sc hlssinger , 1 9 76) uses bulk 
transmission functions for discrete pressure layers. Net 
radiative flux at ground level is a function of incoming 
solar and longwave radiation and surface albedo. 
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NCGAPS 
lar and vec 
model has a 
forecast da 
increase in 
ing eke dia 
the entire 



is run on the CYBER 205 which has high speed sea 
tor pipeline processors. Presently, tne forscas 
wall clock execution time of thirty minutes pe 
y cn a 2.4 X 3.0 degree grid. Any significan 
rhe model’s performance would require redesign 
hatic processes to include more vector! zation fo 
me d er . 
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III. SPECTRA I DECOMPOSITION ROUT INE 

The complex appearance of a typical 500 mb heigh- wave 
pattern often conceals the sources of forecast error. The 
pattern's complexity is due to the superposition of plane- 
tary, long, medium and synoptic waves. Every pattern, 
regardless of its intricacy, can be represented as the sum 
of a series of pure harmonic waves. A spectral decomposi- 
tion routine will resclve a complex wave into its 

components, that is, the amplitude and phase angle of each 
wave number 1, 2, 3, etc. The spectral analysis used in 

this stud j was develojed at the Naval Postgraduate School by 
Prof. R. T. williams. 

The 50C at wave pattern along 42N (Fig. 2) for 18 Febru- 
ary will be spectrally decomposed tc illustrate the 

t echr.ique . 

To facilitate this illustration and the foregoing 
sequence of illustrations. Fig. 3 has an abcissa/or din at e 
scale ratic of 25000k m/60 Cm or approximately 42000 to 1. 

Note the negative longitudes are west and the positive 
lonaitudes are east. A careful examination of the heights 
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500MB SFC HEIGHT 
FEB 18 




Figure 2. 



The Heicrht contours of the 5O0ab Surface a* QOZ 
13 Fet. Contour Intervals are 120 Meters. 



HEIGHT Of 500MB SFC FOR 42N 




LONG ITUOEt DEGREES 3 



Figure 3. The Heigh* of the 500 nib Surface Along 42N at 00 
13 Fee. The Negative and Positive Longitudes a 
West and East Fesp ect iv=-ly of the Greenwich 
Meridian. Boxes Appear at 30 Degree Intervals. 
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along 42N suggest a relatively significant amplitude of 
wavenumber 3 as evidenced by the prominent ridge locations 
in the vicinity of 12QW, 0 and 90S. Also the preponderance 
of ridges and troughs in the 4000km (50-60 degrees) range 

are the signatures of wavenumbers 6 and 7. Since a minimum 
of six points is required to adequately resolve a wave, the 
maximum wavenumber accurately rssclveabls with a 72-point 
series (five degrees of longitude) is 12. A horizontal res- 
olution of three degrees of longitude can accurately depict 
wavenumbers up to 20. Amplitude estimates are difficult and 
finer resolution would be speculative at best. 

1 he numerical results: of the decomposition routine are 
presented in Table 1 ( although only wavenumbers 1-20 can 

be accurately resolved, the tables and figures in this study 
will represent wavenumbers 1 - 28) . The phase angle and 

amplitude (ir. meters about the mean height) of the 
wavenumber 1 through 29 components are shown. 

first the mean height was calculated, then the amplitude 
and phase angle of each wavenumber 1 - 28. The amplitudes 
are giver, as difference values about the mean height and the 
phase angles are in reference to the center of the longitu- 
dinal domain (0 degrees) . 
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TAB L2 I 



The numerical results 

Spectral decomposition of 
meters and the phases are 



cf the spectral decomposition. 

500mb wave at 42N. Amplitudes a 
in degrees. 



ARITHMETIC MEAN 
5485.91 



WAVE NUM 


AMPLITUDE 


PHASE 


1 


52.04 


-88.94 


2 


66.67 


27.38 


3 


92. 1 7 


-157.89 


4 


24.91 


88.19 


5 


26.40 


126.31 


6 


46. 92 


-41.79 


7 


57. 1 7 


129.36 


8 


13. 86 


134.82 


c 


13.63 


-0.63 


10 


10.19 


- 115.47 


1 1 


21.45 


101.96 


12 


10.36 


135.11 


13 


4.26 


-41.04 


14 


1.20 


142.10 


15 


6.06 


111.40 


16 


4.96 


125.36 


17 


9.54 


77.27 


18 


4.70 


1.36 


IS 


4.54 


64.28 


20 


2.9 1 


-52.99 


21 


1.13 


-96.49 


22 


1.56 


27.35 


23 


0.46 


-14.74 


24 


1.30 


- 146.30 


25 


0.4 2 


- 114.46 


26 


1.17 


136.61 


27 


1.30 


41.12 


28 


0.6 5 


-28.40 



A graph of wave number versus amplitude (Fig. 4) f u it 
clarifies the high relative amplitudes of wavenumbers 3, 
and 7 proposed by our initial visual analysis. The ccmpa 
tive importance cf wavenumber 11 is also apparent. 
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SPECTRAL DECOMPOSITION OF 500M8 SFC 




Lours 4. The Relative Anplitudss of the Wavenumber 1-28 

Components of tne SOOmb Surface Aiona 42N at 00Z 
13 Feb. Boxes Appear for each Wavenumber 1-23. 



HEIGHT OF 500MB SFC ANO UNI FOR 42N 




LONGITUDE I DEGREES) 



gure 5. The Height of the oOOmb Surface and Wavenumber 
Components Along 4 2N at 002 on 13 Feb. The 
Nscative and Positive Longitudes are West and 
East Respectively of the Greenwich Meridian. 
Bcxes are Positioned at 30 Degree Intervals, a) 
Wavenumber 1 Component. 
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HEIGHT OF 500MB 5FC AND WN2 FOR 42N 





HEIGHT OF 5 COM3 SFC ANO HN1-3 FOR 42N 
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The series cf graphs in Figs. 5a-e elucidate the signif- 
icance of the plane nary component cf a he 5 00 mb wave. All c£ 
these figures craph the height versus the longitude cf -.he 
original wavs (solid line) with various wavenumber combina- 
tions (dashed line) superposed for comparison. Figs. 5a, b 
and c chart wavenumbers 1, 2 and 3 respectively, while d is 
the summation cf wavenumbers 1, 2 and 3. 



t 




5e) Summation of wavenumber 1-28 components 



Note hew these three waves begin to resemble the origi- 
nal pattern. It is important to realize the summation of 
the full series cf harmonic waves would reconstitute the 500 
mb wave excepting any high frequency variations (greater 
than wavenumber 23) present in the original wave. Frg. 5s 
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500MB WN3 POST FOR 42N NOGRP 0-72 




i i i ^ "”* 7 ” ,i»ii 

-ieo.o-:so.3 -:2C.o -eo.o -so.o -30.c o.o 30.0 so.o so.o 120.0 isa.o 1 a 



LCNS I TUDE l DEGREES ) 



500MB WN3 RNRLY5I5 FOR 42N NOGRP 0-72 




Figure 6. 



Ho vrc 1 l c ~ Diagranis of tha Wavenumber 3 Componer.' 
Sftai500mb Surface Alone 42N. Forecasts aaa 
Analyses ai3 fcr 00Z cn 18-21 Fab irom n r S . 
intervals are 3C Meters. 
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represenes the sum of wavenumbers 1 through 28 and, on this 
scale, appears tc coincide with the test wavs. 

This analysis approach has proved to be a useful 
research tccl when applied to forecast verification (Baum- 
hefner and Downey, 1978; Somerville, 'i 980) . Spectral decompc- 
siticr. cf a forecast anc its verification analysis allow 
direct comparison of the a id lit ude and phase speed of each 
wavenumber. The dimension of time further strengthens this 
analysis technique via a Hovmoll-er diagram. These trough 
and ridge lcngituda-t im e plots along a fixed latitude depict 
-he time svcluticn cf the forecast variables. Fig. 6 com- 
pares the amplitude and phase speed of successive 24, 43 3r.d 
72 hour forecasts for wavenumber 3 with the observed atmos- 
pheric behavicr. 

In the next chapter this and other cases will be ana- 
lyzed in this manner. The harmonics of each forecast are 
grouped ir.tc planets ry (1 - 3) , long (5-7) and medium (8-12) 
waves before comparison to the verification analysis via 
Hovmcller diagrams. Readily discernable errors in ampli- 
tude, phase speed, development of barociinic systems, damp- 
ing cf smaller scale features and planetary wave structure 
are jest a few of the advantages cf verification using spec- 
tral components. 
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IV 



• SPECTRAL VE RIF ICATION CASE STUDIES 

In this chapter the spectral verification presented in 
chapter 3 will he applied to three NCGAPS forecasts, one for 
three days and two extending for five days. All forecasts 
are made from 00Z data anc are verified at 24, 48, 72, 96 

and 120 hour marks of the forecast. 

4.1 Case 1 - February 18 to 21, 1982: 

This period was selected on the bas:.s of its pronounced 
planetary and long wave activi-y. This prediction was made 
with an earlier 'coarse resolution* (4 degrees iat. by 5 
degrees long.) version of the NOGAPS forecast model. The 
spectral decomposition was applied tc features along 42N 
averaged between 38 N and 46N. 

The initial 500mb pattern for case 1 was discussed in 
Chapter 3 and is shown in Fig. 2. Prominent features along 

42N include planetary wave troughs at 160 w and 45 W and 
strong ridges at 170E, 1 10W and 152, suggesting a strong 

wavenumber 3 component amplitude. The preponderance cf 
ridges and troughs with wavelengths cn the order of 50 tc 60 
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degrees are the signa xures of wavenumbers 6 and 7. Recall 
“hat spectral decoaiposit ion of the heights along 42N con- 
firmed the relative dominance cf wavenumber 3, 7 and 11 

components (Fig. 4) . 

The 72-hcur forecast and analysis for the 500mb surface 
on 21 Feb (Figs. 7a and b) appear quite similar under a cur- 
sory inspection. 

The trough at 60-70W has been under forecast by the model, 
but overall the forecast seems to be reasonably successful. 
Note the excellent forecast of the two troughs in the 
Pacific Ocean (180W and 140W) and the strong ridge over 
Europe. 

fchen the spectral decomposition is applied to the 72 
hour forecast and its verification (Fig. 7c) significant 
amplitude differences are present. Specifically, in the 
planetary wave group (wavenumbers 1-3) the observed atmos- 
phere has maintained, although to a lesser degree, the 
dominance cf wavenumber 3 whereas the model has shifted the 
peak amplitude away from wavenumber 3 and into wavenumber 1. 
In the long wave group (wavenumbers 5,6 7) the dominance cf 
wavenumbers 6 and 7 has shifted to the sole peak in 
wavenumber 7 in both the ncdel and the analysis. 
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a > 500MB SFC HEIGHT 

NOGRPS 72 HR TOR FEB 21 




500MB SFC HEIGHT 
FEB 21 




Figure 7. 



50 Gmt Heights and Wavenumber Comparer. 
21 Fst. a) Forecast and o) Anaxysis 
Ir.t ervals cf 120 Meters, c) Wavenumber 
Amplitudes in Meters. 



ts to 
with C 
CCiUDO 



0 



31 



I o o 
C2 M 




7c) Wavenumber component amplitudes. 



NOGAtS accurately forecast the wavenumber amplitudes ir. this 
group. The cedium wave group (wavenumber 8-12 components) 
initially indicates a preference for wavenumber 11, but 
after 72 hours the significance has shifted to the lower 
wavenumbers 8 arc 9. This shift is slightly more pronounced 
in the model resulting in erroneously high rorecas- ampli- 



tudes for wavenumbers 8 aid 9. 

The Hovmoller dragrams for the planetary wave components 
of the forecast and the analysis (Figs. 8a and b) clearly 
illustrate these differences. 
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T I ME ( HOURS ) cr , T ] ME l HOURS ) 



a) 



PW FCST FOR 38-46N NOGRP 0~72 
FEB 18 - FEB 21 




PW RNRLYSIS FOR 38~46N 
FEB 18 - FEB 21 




figure 8. Planetary 'Java Hcvmoller Diagrams for 18-21 ?eb. 
a) rcraca st and b) Analysis wrth Contour 
Intervals of 5C Me-ars/ 



The model has accurately positioned ridges at 120W and 0, 



but 



^3 3 



u n d e r e st 



mated 



the amplitude 



of the former by 50 
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meters and overestimated the latter by 50 meters. The -wc 
planetary wave troughs originating at 180W and 40W are weak- 
ened too rapidly in the forecast and become distorted by the 
spurious formation of a deep trough at 9CE. The trough 
originating at 30E also weakened toe rapidly to the point 
where it is higher than the mean height and is only a trough 
in the sense that it is between two ridges. The ether prom- 
inent feature is the forecast trough at TOE which has nc 
counterpart in the observations. The seemingly random 
changes in trough an c ridge intensity form a pattern when 
the differences in component amplitudes are considered. 

The NCGAFS 72-hour forecast shews the wavenumber 1 sig- 
nature emerging with general ridging from 150W to 40E and 
general troughir.g from 40E to 1 50 v? while the analysis 
depicts mere of a wavenumber 2 and 3 signature. Figs. 9a 
and b compare the forecast and observed wavenumbers 1 and 3 
ir.di vidua lly. 

NOGAFS forecasts wavenumber 1 to be 70S more intense than 
the observation while under forecast ing wavenumber 3 by 55S. 
Despite the presence of wavenumbers 2 and 3 , the planetary 
wave char ec* eristics in the 72 hour NOGAPS forecast are 
largely determined by the overwhelming dominance of 
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ANALYSIS UN3 COMPONENT ISULiUJ 
fCRECfiST WN3 COMPONENT I08SHED) 




Figure 9. Wavenumber 1 and 3 Components for Cass 1. a) 
Wavsr.umber 1 Ccirponent for Analysis (Solid) an 
Forecast (Cashed) for 003 21 Fso. b) Sams as 
ex cep- for Wavs number 3 Co aeons nr. 
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strong, exist in the regicn of 120W to 6QE. Synoptic exam- 
ples of the effect of the reduction of the wavenumber 3 
component are evident in the ridge at 120 W and the trough at 
SOW. From 150W to 1 10W the model under forecasts the ridge 
in this area and smooths cut its amplitude. From 110W to 
SOW the model over f o reca sts the height of the eastern side 
of the ridge and smooths out the trough at SOW. These 
amplitude errors are also evident in other regions of the 72 
hour forecast (Figs. 7a and b) . 

The long wave group verifies the wavenumber 5,6 and 7 
components. The general agreement between the individual 
wavenumber amplitudes in this group (Fig. 7c) through the 
72-hcur period consequently produces a general agreement in 
the amplitudes of their respective Hovmolier diagrams (Figs. 
10a and b ) . 

This diagram shews seven troughs ar. d ridges at 72 hours sig- 
nifying the dominance of wavenumber 7. in this grouping a 
regicn of maximum amplitude (contour intervals are 15m) is 
evident between 180W and 30W in beta the model and the 
observation. There is strong agreement between the forecast 
and observed long waves in the Hovmolier diagrams shewing 
NOGAFS performing well at this scale. 
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a) LW FCST FOR 38-46N NQGflP 0-72 

' FEB 18 - FEB 21 




b) LW RNRLYSIS FOR 38-46N 

FEB 18 - FEB 21 




Figure 10. Iona Wave Hovmcllsr Diagrams for 18-21 Feb. a) 
Forecast and t) Analysis with Contour Interval* 5 
cf 30 Meters. 



The interaction of wavenum ber s 8 through 12 in the 
medium wave group creates a complex pattern in the analysis 
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and forecast trough-ridge diagrams. The Hovmoiier diagrams 
for the forecast and observation (Figs. 1 la ana b) both show 
eleven waves initially and nine waves after 72 hours. 

The amplitudes (contour intervals are 10 meters) ar. d phase 
speeds can still be compared, but the errors are more diffi- 
cult to trace. The model shows an erroneous loss of ampli- 
tude, particularly evident at 48 hours. However, at 72 
hours the increased activity in the region from 140E to 12QW 
is properly forecast. The phase speed is approximately 10 
degrees tcc fast in the region from 00E to 120E. Con- 
versely, the model has done extremely well in forecasting 
the speed cf the wave features between 14GE and I50W, Over- 
all, the model exhibits reasonable success for both the 
amplitudes and phase speeds cf the medium waves at 7 2 hours. 
It is important to note the amplitude cf the medium wavenum- 
ber s is often smaller than the planetary wavenumbers by as 
much as an order of magnitude. Also as wavenumbers increase 
to the point where they cannot be accurately represented by 
the model's horizontal resolution, the spectral decomposi- 
tion is less accurate as well. 

The results of Case 1 indicate NOGAPS erroneously 
shifted the planetary wavenumber dominance from wavs 3 and 
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TIME (HOURS) cr TIME(HOURS) 



a) 



MW PCST FOR 38-46N NOGRP 0-72 
FEB 18 - FEB 21 




m RNRLYSIS FOR 38-46N 
FEB 18 - FEB 21' 




Figure 11. Medium Wavs Hcvmcller Diagrams for 18-21 Fsb. 

a) Fc recast aid b) Analysis with Contour 
Intervals of 15 Meters. 
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model and the observation were difficult. Both the ampli- 
tudes and phase speeds of the long waves were accurately 
predicted by the model. NOGAPS shewed some less cf ampli- 
tude in the medium waves and shifted peak amplitude from 
wavenumber 11 to 8 and 9. The phase speed errors of the 
medium waves as a group were small. 

4.2 Case 2 - October 06 tc 11, 1982: 

This period provides a good example of the formation cf 
an intense ridge (blocking high) over the western U.s. and 
was discussed in the NOAA, western Regional Attachment of 
October, 1 962. For this case NOG APS horizontal resolution is 
2.4 degrees of latitude by 3. degrees of longitude interpo- 
lated onto a 2.5 by 2.5 degree so her real grid. The spectral 
decompcsitic r. is applied tc the height of the 500mb surface 
along 45N (averaged between the heights given at 42. 5N and 
47 . 5 N ) . 

The 50Cmb heights, depicted from 30N to 60N, for 06 Cct 
in Fig. 12a shows a fairly zonal pattern from the east coast 
of Asia tc the Eastern Atlantic. 

Note the ordinate/abscissa ratio is more than one. This 
distortion produces a visual magnification cf the amplitudes 
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a) 500MB SFC HEIGHT 

OCT 06 




SPECTRAL DECOMPOSITION OF 500M8 SFC 
ANALYSIS 06 OCT AT 45N 




Figure 12. 



5CCmt Heiahts ana Wavenumber Components for 06 
Cot. a) 500m t Heights with Contour Intervals of 
120 Meters, b ) Wavenu mber Component Amplitudes in 
Meters. 
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a-: 10W. Spectral decomposition (Fig. 12b) shews wavenumber 
2, 5 and 7 components are dominant initially. The zonal 

character cf this case is shown by the relatively even dis- 
tribution cf amplitudes ancng wavenumbers 1 through 8. 

Seventy-two hours later (Fig. 13a) a deep trough has 
developed ever the central U.S. and prominent ridges have 
rapidly developed over the western and eastern U.S. coasts. 
The model (Fig. 13b) dees not adequately forecast the 
strength cf these features. A comparison of the spectral 
decomposition (Fig. 13c) cf the analysts and the forecast at 
72 heurs shows a marked increase in the dominance cf the 
wavenumber 2 anc 7 ccmponents cf the atmosphere while the 
model incorrectly smoothed the amplitude of wavenumber 2 but 
correctly maintained the relative strength of wavenumber 7. 

At the 120 hour mark (11 Oct) of the forecast period 
both the eastern and western U.S. ridges have continued tc 
strengthen (Fic. 14a) . The model shows an ability tc prop- 
erly forecast the phase speed of the wave features such as 
the troughs at 170W and 40W and the ridge at 160S but con- 
tinues to underestimate the strength of the ridges (Fig. 
1 4 b ) over the U.S. Spectral decomposition of the observa- 
tion shows the largest amplitudes are in wavenumbers 2, 3, 4 
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a) 500MB SFC HEIGHT 

NOGRPS 72 HR FOR OCT 09 

60* N 
55* N 
50* M 
«* N 
40* N 
35* N 
30* N 

I 

b) 500MB SFC HEIGHT 

OCT 09 

60* N 
55* N 
SO* N 
15* X 
10 * H 
35* N 
30* N- 

60* E SO* £ 120* E 150* E 180* M ISO* » 120* H 90* M 60* M 30* H O' 30* E « 



SPECTRAL DECOMPOSITION OP 500MB SFC 
RNRLYSISISOLID) AND 72HR FORECAST ( DfiSHED ) 
09 OCT AT i5H 
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and 5 all having similar valuas (Fig. 14c). Meanwhile the 
model has shifted the dominance to wavenumber 4 while keep- 
ing the amplitude of wavenumbers 5 and greater relatively 
small. The largest amplitude changes in the model during 
the forecast period occurred in wavenumbers 2, 3 and 4 while 
the atmosphere additionally showed large amplitude changes 
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planetary wave trough-ridge diagrams (Figs. 15a and 



exhibit the 
cc mpcnen t . 


characteristics of a large wavenumber 2 


After / 2 


hours the observation has greatly increased the 


strength of 


the wavenumber 2 signature while the forecast 


decreased it 


(Fig. 16a). The model has successfully fcre- 


cast the st: 


rengtfc of the trough at 1S0W (-135m), but the 


strength cf 


the other trough and both ridges are severely 



underestimated. The combined dominance of wavenumbers 2 and 
3 in the NCGAES forecast has introduced a spurious ridge at 
OW and trough at 50E. 

Ey 120 hours the atmosphere exhibits maximum amplitude 
in the planetary wavenumbers 2 and 3 (Fig. 14c). This shift 
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500MB 5FC HEIGHT 
NOGRPS 120 HR FOR OCT 11 




b) 500MB 5FC HEIGHT 

OCT 11 




SPECTRAL DECOMPOSITION OF 500MB SFC 
ANALYSIS (SOLID) AND 120HR FORECAST ( DASHED ) 




Figure 14. Same as Fig. 7 except for 11 Oct. 
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wavenumber 3 signature. The model has forecast the ampli- 
tude of wavenumber 1 accurately while underestimating 
wavenumber 3 by 30% (Fig. 16b). 

Thus tbs model has ones again failed to predict the 
shift and erroneously smoothed the highest amplitudes. The 
model's planetary wave Hcvmoller diagram shows the weaker 
amplitudes at 120 hours and, significantly, the resultant 
smoother pattern. Specifically the ridges at 105E and 110W 
are too lew and the trough at 45W is toe high. Synopticsliy 
these errors are evident in Fig. 14b as the failure of the 
model to adequately forecast the strength of the ridge 
building on the west coast of the U.S. and the trough at 
40 W. The general placement of comparable features is good 
especially in the Pacific region near the dateline. 

The long wave component Hcvmoller diagrams (Figs.' 17a 
and t) shew NCGAPS successfully forecasting the long wave 
features. Overall the forecast has accurately predicted the 
general regions of high long wave amplitudes, 120 V to 90E, 

and lew long wave amplitudes, at 90S to 1204. observed long 
wave amplitudes are more intense than the forecast after the 
72 heur mark. Comparison cf model and observed phase speeds 
shows NOGAES to be fast, particularly in the high amplitude 
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HEIGHT (METERS I w HEIGHT(MCTERS) — HEIGHT l METERS) 

5200.0 S10Q.0 5600.0 5000.0 5200.0 5100.0 5600.0 5600.0 5200.0 5100.0 5600.0 5800.0 



ANALYSIS WN2 COMPONENT ISOLlDT 
FORECAST HN2 COMPONENT (DASHED) 
72 HOURS 109 OCT) 




“ANALYSIS WN3 ' C OHPOnENTTSOCI D )’ 
FORECAST WN3 COMPONENT (DASHED) 
120 HOURS 111 OCT) 




ANALYSIS WN5 COMPONENT (SOLID) 
FORECAST WN5 COMPONENT (OASHEO) 
120 HOURS 111 OCT) 




Figure 16. Wavenumber 2, 3 a r.d 5 Components for Case 2. 

Same as Fig. 9 except a) Wavenumber 2 for 09 
Oct, b) Wavenumber 3 for 11 Oct and c) 
Wavenumber 5 for 11 Oct. 
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LW FCST AT 42.5 - 47. 5N 0-120HR 
OCT 06 - OCT 11 




L U ANALYSIS AT 42.5 - 47. 5N 
OCT 06 - OCT 11 
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MW FCST FOR 42.5 - 47. 5N 0-1 20H 
OCT 06 - OCT 11 




MW ANALYSIS AT 42.5 - 47. 5N 
OCT 06 - OCT 11 
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direct comparisons car be made. Sy r.opt icaily this differ- 
ence is net visible cn the scale of Figs. 14a and b but 
would shew up as a slightly rapid movement of features with 
a wavelength cf 30 to 45 degrees over the 120 hour forecast 
per icd. 

The results of Case 2 indicate NOGAPS erroneously 
smoothed a shift in the planetary wave activity at 72 hours 
and again at 120 hours. Planetary wave phase speed errors 
were small fer comparable features and negligible in the 
Pacific region. Although long wave forecast amplitudes were 
smoothed, NCGAPS correctly produced twe distinct regions cf 
contrasting amplitudes. Ehase speeds were slightly too fast 
in the high amplitude region but satisfactory in the low 
amplitude regicn. The amplitudes cf the medium wave fea- 
tures were well forecast, but shewed a fast bias in wave 
speed and several hey features were lest in the forecast. 

4.3 Case 3 - Cctcber 22 tc 27, 1982: 

Curing this forecast period the 500mb surface changed 
frem a predominantly long wave flow pattern to a more zonal 
situation. In particular the region from 150F to 140W 
underwent the breakdown cf a blocking ridge. Again the 
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NOGAES horizontal resolution is 2.4 degrees la-. by 3 
degrees long. interpolated onto a 2.5 by 2.5 degree spheri- 
cal grid. The spectral decoin position is applied to the 
heights alone 45N averaged between 42. 5N and 47. 5N. 

The 500mb heights from 30N to 60N for 22 Oct (Fig. 19a) 
includes prominent troughs at 6QE, 125E and 14QW with lesser 
troughs at 7 5 and 1 5 W . 

Noteahie ridges exist at 90E, 170E, 110W, 45W and 20E. The 
wavenumber 5 dominance is clearly present in Fig. 19b alone 
with strong amplitude components of wavenumbers 3 and 4 due 
to their particularly strong trough-ridge systems. 

After 72 hours amplitude differences are already appar- 
ent between the 500mb patterns of the model and the analysis 
(Figs. 20a and b) . The model has underforecast the strength 
of the ridges at 90E and 35E and the troughs at 70S, 140E, 
1 4 0 W and ICE. The only obvious NOGAES overforecast is in 
the ridge at 100W. 

The individual wavenumber components (Fig. 20c) a- 72 
hours show the dominance remains with 4 and 5. Large ampli- 
tude increases have occurred in wavenumbers 1, 7, 9 and 13 
and large decreases at 2 and 3. The model has accurately 
predicted most of the changes but consistently 
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b) 



SPECTRAL DECOMPOSITION OF 500MB SFC 
ANALYSIS 22 OCT AT 45N 




underestimates the strength of wavenumbers 6 through 10. In 
terms of this study's wave groupings the amplitudes of the 
planetary waves are forecast very well, the long waves are 
too weak ty 10 to 203 and the medium waves are too weak fcv 
30 tc U03. 
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500MB SFC HEIGHT 
NOGRPS 72 HR FOR OCT 25 




b) 



500MB SFC HEIGHT 




SPECTRfiL DECOMPOSITION OF SOOMB SFC 
ANALYSIS ( SOLID J AND 72HR FORECAST IDRSHED) 
25 OCT RT 4SN 
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At the 120 hour mark both pha se speed ar.d amplitude 
errors are evident in the 500mb patterns (Figs. 2 la ar.d fc) . 
Phase speeds are consists r.tly tcc rapid. Examples include 
the trough at 75E (10 degrees of longitude too fast), the 
ridge at 180W (15 degrees too fast) and the trough at aOW 
(25 degrees tcc fast). The amplitudes of the trough at 160E 
and the ridge at 1 8 0 W are too weak and the ridge at 95 W is 
toe strong. A ridge has formed in the region from 160W to 
120W which has aided in the weakening of the blocking ridge 
at 160W. The model has missed this important feature 
en ti rel y . 

The individual wavenumber components (Fig. 21c) at the 
120 hour mark show a general decrease in all amplitudes with 
especially large percentage drops in wavenumbers 3, 4, 5, 6 

and 1. The relative dominance of planetary waves is much 
smaller new, hut there exists a distinct peak at wavenumber 
4. The model's wavenumber components also show a general 
decrease in amplitudes with significant errors occurring in 
3, 4, 8, 9 and 12 while the dominant wavenumber is 

wavenumber 3. In terms of this study's wave groupings the 
amplitudes of the planetary waves are too strong by 20 to 
30%, the long waves are well forecast and the medium waves 
are too weak by 30 to 40%. 
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a) 



500MB SFC HEIGHT 
NOGRPS 120 HR FOR OCT 27 




b) 500MB SFC HEIGHT 

OCT 27 




SPECTRAL DECOMPOSITION OF 500MB SFC 
ANALYSIS* SOLID 1 AND 120HR FORECAST (CASHED] 
27 OCT AT 45N 




Figure 21. Same as Fig. 7 except for 27 Occ 
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The Hcvmciler diagrams for ths planetary wave components 
of ths forecast and the analysis (Figs. 22a and b) both show 
a shift from the sole dominance of wavenumber 3 to a shared 
dominance with wavenumber 1 by 72 hours. 

The forecast phase speed error of the model is negligible at 
this point. The wavenumber 3 features are forecast too 
strong as a result cf the 20% ovarforecast by the model. 
Examples include the troughs at 1102 and 150W and the ridges 
at 80S and 30E. At 120 hours the model shifts its maximum 
activity tc wavenumber 3 while the atmosphere shows a dis- 
tinct wavenumber 1 signature. Now the NOGAPS wavenumber 3 
features, the troughs at 1202, 50W and 35E and the ridges at 
85W and 35E, are much toe strong. Fig. 23 graphically pres- 
ents the 100% eve rf or sc as t of the wavenumber 3 component by 
the model. 

The positions of the ridges at 85W and 352 are excellent, 
but the trough at 1202 has been erroneously retrograded 15 
degrees. Another interesting aspect of this case is the 
retrogression cf the planetary wave which originally was at 
150W. NOGAPS handles this phase shift very well. 

Within ths long wave group the signature of wavenumber 5 
persists through 72 hours in both the model and analysis 
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a) 



PW PCST AT 42.5 - 47. 5N 0-120HR 
OCT 22 - OCT 27 




PW ANALYSIS AT 42.5 - 47. 5N 
OCT 22 - OCT 27 




(Fiqs. 24a anc t) . A general comparison of observed and 

forecast lcng waves shows good agreement on this scale, 
particularly in phase speeds. The model urderf or ecasts ohe 
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"ANALYSIS WNI TTOiiPONENT l SOL 10) 
FORECAST WN3 COMPONENT 10RSHED) 
120 HOURS 127 OCT) 




Figure 23. Wavenumber 3 component for Casa 3. Same as Fig. 
9 except. Wavenumber 3 for 27 Oct. 



strength cf wavenumbers 6 and 7 producing slightly weak 
amplitudes in the region from 302 to 1202. 

However, by 120 hours both observed and forecast pat- 
terns show a cuel dominance in wavenumbers 5 ana 6 with the 
model amplitude a little too weak in wavenumber 5 and a lit- 
tle too strong in wavenumber 6. In the last 24 hours of the 
forecast pericc the atmosphere introduced a sixth trough and 
ridge by a general retrograding of certain existing features 
(1S0W to 4CE) by 5 tc 1C degrees of longitude. The model 
introduced a sixth trough and ridge by a rapid eastward 
movement cf certain features (25W to 140E) . This incorrect 
method cf handling the dominance shift resulted in phase 
speed errors in the range of 5 to 15 degrees too fast. 
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TIME (HOURS) ^ TIME (HOURS) 



a) LH ECST RT 42.5 - 47. 5N 0-120HR 

OCT 22 - OCT 27 




LW ANALYSIS RT 42.5 - 47. 5N 
OCT 22 - OCT 27 




In the medium wave group (Figs. 25a and b) the model 
correctly began with the signature of wavenumber 9, main- 
rained i t through 72 hours then shifted to wavenumber 8 
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a) MW FCST FOR 42.5 - 47. 5N 0-120H 

OCT 22 - OCT 27 




b) MW RNfOSIS AT 42.5 - 47. 5N 

OCT 22 - OCT 27 




dominance ty 120 hours. 3cth diagrams show rheir largesr 
amplitudes in the region from 80 E no 120W. Overall rhe mod- 
el's amplitudes are too weak as a resuir of consist enrly 
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underfcrecasrirg the strength of the medium waves. The 
phase speeds are too fast by apprcxi 
region from 2CE tc 16 CE. 
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• CCNC1USICNS J| SUGGESTIONS FOR FUTURE 

In this study one three-day and two five-day NOGAPS 500®t 
forecasts were spectrally verified. The cases were chosen 
to reflect three different planetary ana long wave situ- 
ations. In Case 1 the SCCmb pattern is characterized by a 
via c reus three wave pattern. In Case 2 a formation of a 
blocking high and a shift to a low index flow dominate the 
500ah pattern while Case 3 describes the breakdown of a 
block and a shift to a high index flow. 

Two analysis techniques were used to study the observed 
and forecast spectral statistics. The amplitudes of each 
wavenumber in the analysis and the forecast are compared at 
key periods during the forecast and Hovmoller diagrams are 
prepared fer key wavenumber groupings following Baumhsfner 
and Ecvney (1978). Save number components were grouped into 

planetary (waves 1-3) long (waves 5-7) and medium waves 
(waves 8-12) and displayed or. trough-ridge diagrams. 
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failed io correctly predict large amplitude charges ir the 
individual planetary wavenumber components. NOGAPS both 
predicted large amplitude changes which did not occur in the 
observations and failed to predict large amplitude changes 
which dad occur in the analysis. In Case 1 the model fore- 
cast an erroneously large amplitude increase in wavenumber 1 
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The medium waves (6-12) verification exhibited some sys- 
tematic errors. HOG APS tended to smooth the amplitudes cf 

the medium waves to the point where key features were often 
lost. Also phase speeds cf the medium waves were generally 
too fast in the model. Specifically, examples of fast phase 
speeds in Case 1 are in the region from OE to 129S, in Case 
2 from 1 9 0 5i tc 30W and in Case 3 from 20E to 160E. 

Several suggestions for future research follow from this 
study. First, the number of forecasts spectrally verified 
should be increased. New that the spectral verification 
software has been prepared, the verification could be done 
routinely. 1 he coal of the increased number of cases would 

be tc document systematic spectral errors in NOGAPS. 

A second area that should be investigated is ■‘•he rela- 
tionship between *:he strength cf the circumpolar vortex and 
the forecast phase speed errors. Variations in the 
wavenumber zero component can be found by comparing gradi- 



ents cf the mean 


height of the 


forecast 


tc z he coserva 


The strength cf 


the circumpolar 


vortex 


may be a taez 


the phase speeds 


of the shorter 


synoptic 


waves. 
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A third suggestion involves the question cf the optimum 
north-south interval for verification. This study used 
eight degrees of latitude for Case 1 and five degrees for 
Cases 2 and 3. Further studies should examine the effects 
of narrower or wider latitude bands in the analysis. Simi- 
larly the effects of latitude variations on wavenumbers 
could be determined by using data from different latitudes 
(6 ON , 3 ON , etc) . 



This study illustrated the utility of diagnostic verifi- 
cation techniques for study of operational numerical pre- 
diction medals. Future improvement in NOGAPS forecast 
skill, particularly in the medium range, is dependent in 



part cr the continuing study of the model to isolate 
strength and weaknesses. 
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